Water motion and pH jointly impact the availability of dissolved inorganic carbon to macroalgae

The supply of dissolved inorganic carbon to seaweeds is a key factor regulating photosynthesis. Thinner diffusive boundary layers at the seaweed surface or greater seawater carbon dioxide (CO2) concentrations increase CO2 supply to the seaweed surface. This may benefit seaweeds by alleviating carbon limitation either via an increased supply of CO2 that is taken up by passive diffusion, or via the down-regulation of active carbon concentrating mechanisms (CCMs) that enable the utilization of the abundant ion bicarbonate (HCO3−). Laboratory experiments showed that a 5 times increase in water motion increases DIC uptake efficiency in both a non-CCM (Hymenena palmata, Rhodophyta) and CCM (Xiphophora gladiata, Phaeophyceae) seaweed. In a field survey, brown and green seaweeds with active-CCMs maintained their CCM activity under diverse conditions of water motion. Whereas red seaweeds exhibited flexible photosynthetic rates depending on CO2 availability, and species switched from a non-CCM strategy in wave-exposed sites to an active-CCM strategy in sheltered sites where mass transfer of CO2 would be reduced. 97–99% of the seaweed assemblages at both wave-sheltered and exposed sites consisted of active-CCM species. Variable sensitivities to external CO2 would drive different responses to increasing CO2 availability, although dominance of the CCM-strategy suggests this will have minimal impact within shallow seaweed assemblages.


Results
Site information. Gypsum blocks were used to compare the rates of mass transfer at the wave-exposed and wave-sheltered sites. The average mass transfer rates were over two-times higher at the wave-exposed field sites compared to at the wave-sheltered sites (Table 1). Light attenuation, temperature, salinity and nitrate concentrations did not significantly differ between the wave-sheltered and -exposed sites, although there was slightly higher concentration of phosphate and ammonium at the wave-exposed sites (Table 1).
Prevalence of carbon-use strategies in sub-canopy assemblages. Carbon stable isotope signatures (δ 13 C) were used to indicate whether macroalgae had an active CCM are were using HCO 3 − (δ 13 C ≥ − 29‰) or were not using a CCM and relied solely on diffusive CO2 uptake (δ 13 C < − 29‰) 11,20,34,40,41 .The abundance of non-CCM species was less than 3.7% (SE = 0.99, n = 3) of the macroalgal cover within wave-exposed assemblages (Fig. 1b), however, this was significantly more than the wave-sheltered sites where there was only one non-CCM Table 1. Abiotic site measurements of the three wave-exposed and four wave-sheltered sites. Values represent mean ± 95% confidence interval.
Site measurements Wave exposed Wave sheltered  Figure 1. (a) Percent cover of functional groups at wave-exposed (grey) and wave-sheltered (white) sites in macroalgal assemblages. Bars represent mean ± 1 SE (n (exposed) = 3, n (sheltered) = 4). (b) Percent cover of CCM, non-CCM and coralline species at the wave exposed and wave sheltered assemblages. (c) 13 C signatures (δ 13 C) of species present at the two water motion regimes. Species are grouped into functional groups. δ 13 C values that fall below the dashed line at − 29‰ are considered species that solely use carbon dioxide (CO 2 ). (d) 13 C signatures for individual species at the wave-exposed and sheltered sites. Points represent mean ± 1 SE (n = 3 to 21). Asterisks (*) indicate a significant difference between water motion regimes (P < 0.05), tested with a Nested ANOVA (a-c), or one-way ANOVA (d). www.nature.com/scientificreports/ species present at one of the sites (Caulerpa brownii, Chlorophyta), making up 0.05% of the total macroalgal cover at sheltered sites (ANOVA, F 1,5 = 59.3, P < 0.001). The majority of macroalgae possessed δ 13 C value that indicated an active CCM at both the wave-sheltered (mean = 56.3%, SE = 2.33, n = 4) and the wave-exposed sites (mean = 41.49%, SE = 2.83, n = 3; Fig. 1b). Calcareous coralline algae and bare substrate made up the remaining benthic cover (Fig. 1a,b). Three species at the wave-exposed sites, all belonging to the fleshy red functional group, consistently exhibited δ 13 C less than − 29‰ and were thus classified as non-CCM macroalgae: Ectophora depressa (formerly Callophyllis depressa), Hymenena durvillei and Streblocladia glomerulata (Fig. 1c). Both Rhodophyllis membranacea (formerly Rhodophyllis gunnii) and Antithamnionella adnata were fleshy red algae that exhibited varying carbon-use strategies ( Fig. 1d). At the wave-exposed sites, A. adnata and R. membranacea exhibited δ 13 C below − 29‰ indicating sole CO 2 use, however at wave-sheltered sites the δ 13 C of the same species were significantly more positive, up to − 18.8‰ for A. adnata and − 15.8‰ for R. membranacea suggesting active uptake of HCO 3 − (Fig. 1d). The green macroalgae Caulerpa brownii was the only species with a δ 13 C value more negative than − 29‰ at the wave-sheltered sites, and this species was observed at just one of the four wave-sheltered sites.
DIC uptake laboratory experiment. In order to understand how the interactive effects of water motion and an increase in the concentration of CO 2 caused by lowering pH affect CCM and non CCM species, a photosynthesis vs. DIC uptake experiment was conducted under constant light conditions. Water motion, which was simulated by creating a high-or low-mixed environment, and pH both significantly affected the uptake efficiency of DIC (k 0.5(DIC) ) for the non-CCM H. palmata, and the CCM species X. gladiata (Figs. 2, 3). At the present-day pH of 8.1, the k 0.5 of the non-CCM H. palamata was significantly lower in the high-mixed treatment www.nature.com/scientificreports/ compared with the low-mixed treatment (Simultaneous t tests, t = 7.9, P < 0.001), which corresponded to a 25% increase in its DIC uptake efficiency ( Table 2). The effect of mixing was even more significant within the high CO 2 , low pH treatment of 7.6 (Simultaneous t tests, t = 6.7, P < 0.001), with the DIC uptake efficiency of H. palmata being 43% more efficient in high-compared with low-mixed conditions ( Table 2). The lower seawater pH treatment had a concentration of CO 2 that was 230% higher than that of the ambient pH (8.1) treatment ( Fig. 4) and had a similar effect as water motion on increasing the DIC uptake efficiency of H. palmata. The k 0. 5 Table 2. Carbon uptake efficiency (k 0.5 ) and photosynthetic maximum (P max ) of a non-CCM red (Hymenena palmata) and CCM brown species (Xiphophora gladiata) under high-and low-mixing conditions and high (8.1) and low (7.6) pH. Values represent mean ± SE, n = 4. Simultaneous t tests were used to test for significant differences between the water motion and pH treatments within each species, with a p-value < 0.05 indicating a significant difference. Different lowercase letters indicate significant differences between values, while the same letter indicates there is no significant difference.

Species Mixing pH k 0.5 (mM DIC ) P max (μmol O g −1 FW s −1 )
Hymenena palmate (Non-CCM)  Table 2) at a pH of 7.6 under high mixing. In low-mixed conditions the high CO 2 , low pH still resulted in a significantly lower k 0.5 value (Simultaneous t-tests, t = 7.6 P < 0.001), however, the difference was less extreme with a 32% difference between the pH 8.1 and 7.6 treatments ( Table 2). The photosynthetic maximum (P max ) of H. palmata was significantly higher in the high mixed versus the low-mixed treatment at pH 8.1 (Simultaneous t tests, t = 2.88, P = 0.03) but there was no significant difference in the P max at the lower pH treatment of 7.6 even though the CO 2 concentration was 230% higher ( Table 2). The DIC uptake efficiency of X. gladiata also significantly increased in the high mixing treatment (Simultaneous t tests, pH 8.1: t = 18.4, P < 0.001; pH 7.6: t = 40.1, P < 0.001), with a 52 and 69% reduction in the k 0.5 value in the pH 8.1 and 7.6 treatments respectively ( Table 2). The response of the DIC uptake efficiency of X. gladiata to low pH, high CO 2 conditions contrasted, however, to that of H. palmata, with the uptake efficiency of X. gladiata being negatively impacted within the low pH, high CO 2 treatment. Within the high-mixed treatment, the k 0.5 of X. gladiata increased from 0.14 mM to 0.21 mM (Table 2), which corresponded to a significant 50% decrease in uptake efficiency (Simultaneous t tests, t = 8.75, P < 0.001). Within the low-mixed treatment, the negative effect of the low pH treatment was exacerbated even more, showing a 131% decrease in DIC uptake efficiency of X. gladiata (Simultaneous t-tests, t = 33.27, P < 0.001; Table 2). Although the DIC uptake efficiency of X. gladiata was significantly affected by both mixing and pH (Table 2), there was no significant effect to its P max ( Fig. 3; Table 2).
The photosynthesis of both algal species saturated well-below the concentration of DIC in natural seawater of ~ 2 mM: that of X. gladiata saturated at a DIC concentration between 0.3 and 1.36 mM (Fig. 3) and the photosynthesis of H. palmata saturated between 0.41 and 1.09 mM DIC (Fig. 2). A higher variability was observed in the low-mixed compared to the high-mixed treatment for both species. This variability is likely to be a consequence of unstable readings from the O 2 sensor in the low-mixed environment.

Discussion
Water motion regulates the uptake of DIC by macroalgae, with increased mixing enabling more efficient uptake of DIC for photosynthesis. This physiological response was evident in laboratory experiments under controlled conditions, however, the effect on natural macroalgal assemblages in situ was less clear. It was hypothesized that  . The concentrations of CO 2 , CO 3 2− and HCO 3 − in the ambient pH 8.1 and the low pH 7.6 treatment from the DIC uptake laboratory experiments. The low pH conditions cause a shift in the speciation to a higher concentration of CO 2 and a lower concentration CO 3 2− within the seawater. Values were calculated using the 'seacarb' package 42 in R, using a temperature of 12 °C and a salinity of 35 psu, with the assumption that the photosynthesis of the ~ 1 g individuals during the incubations did not have a significant impact on the carbonate chemistry within the chambers. www.nature.com/scientificreports/ a greater availability of CO 2 at wave-exposed sites would lead to more non-CCM macroalgae that rely solely on diffusive CO 2 uptake compared to wave-sheltered sites. A greater abundance and diversity of non-CCM macroalgae was found at the wave-exposed versus the wave-sheltered sites, however, this accounted for only 3 species with a total benthic cover of 3.7%. When we examine this in situ response at the functional group level, however, we observe that fleshy red macroalgae, which are a group where the non-CCM uptake strategy is most common 11 , respond strongly to the water motion regime compared to brown and green macroalgae. The varying responses between functional groups in the laboratory compared to the field suggest that water motion is but one of many factors regulating the uptake of DIC by macroalgae in situ. Alongside water motion, the light regime, nutrient availability, the life strategy of a species and the plasticity of a species carbon-uptake strategy all can influence carbon uptake in natural macroalgal assemblages, making it difficult to predict how entire macroalgal assemblages will respond to changes in the seawater carbon chemistry. The majority of macroalgal species within the kelp forest assemblages of southern New Zealand possessed values of δ 13 C indicative of an active CCM 11,41 at both wave-sheltered and -exposed sites. This dominance suggests that the carbon demands of most of the species in these assemblages cannot be sustained by CO 2 alone and additional DIC from the large HCO 3 − pool is required 5,43 i.e. the photosynthetic rates of the macroalgae exceed what can be sustained with only diffusive uptake of CO 2 . Active CCMs are almost universal in both brown and green macroalgae, with reports of only a handful of green macroalgal species exhibiting δ 13 C less than − 29 indicating sole-CO 2 use 11,34,36,44,45 . These non-CCM brown and green macroalgae are typically found in low light environments, either in deep habitats or beneath dense canopies, where the carbon-demands of individuals are lower 20,34,36 . Given browns followed by greens were the dominant functional groups at both the shallow waveexposed and -sheltered sites, the high abundance of the CCM strategy is expected.
An interesting observation was the difference in within-species variability of the δ 13 C. The dominant canopy forming brown species, M. pyrifera showed very little variability in carbon-use between individuals at different sites, a pattern consistent with previous work 24 . M. pyrifera has an efficient mechanism of bicarbonate uptake, via an anion-exchange protein 10,46 , which allows it to maintain high growth rates 24 irrespective of the availability of CO 2 . Other brown and green species that occurred at both water motion regimes showed more variability in their δ 13 C while still exhibiting HCO 3 − use. Given the sites are constantly subtidal and are all linked by the same body of water with no nearby freshwater inputs, we assume that the 13 C/ 12 C source is the same at all sites 11 . The lower δ 13 C values could suggest downregulation of the CCM, however, 13 C signatures on their own cannot be used to disentangle the variability in 13 C signatures that are more positive than − 29‰ 11 . Local differences in the light environment within the kelp forest canopy, CO 2 availability through water motion and differences in the health and size of individuals and their acclimation to the light environment could all contribute to different carbon requirements and variability of the δ 13 C values from the active-CCM species.
The greatest within-species variability of carbon-use was observed in the fleshy red species. The Rhodophyta phylum has the highest abundance of species with δ 13 C less than − 29, with 35.6% of species that have been analysed showing δ 13 C indicative of sole-CO 2 use 11 . Fleshy red species, Antithamnionella adnata and Rhodophyllis membranacea, were observed in this study to have highly variable carbon-use strategies, both possessing 13 C signatures indicative of HCO 3 − use at the wave-sheltered sites, however, having more negative 13 C signatures at the wave-exposed sites that suggested the sole use of CO 2 . Modulation of CCM activity has been extensively studied within the green unicellular alga Chlamydomonas reinhardtii [47][48][49] and has also been identified in Ulva lactuca 50 . Although modulation of the CCM by algae from the Rhodophyta phylum has not yet been examined, the species we observed seem to possess an ability to optimise their carbon-use depending upon the external CO 2 availability. Such a mechanism would allow these algae to minimise their energy expenditure, but also makes them sensitive to the external CO 2 availability.
By increasing rates of advective transport and reducing the thickness of boundary layers 25,51 , we observe in the laboratory experiments that water motion allows macroalgae to take up DIC at a faster rate. This physiological effect of an increase in the availability of DIC does not appear to strongly impact brown and green macroalgae with active CCMs, as they are likely to be saturated for carbon at current DIC concentrations and instead their photosynthetic rate is more limited by other factors, such as light or nutrients. We observed this DIC saturation in the experiments with X. gladiata, which reached its maximum photosynthetic rate (P max ) at DIC levels below half that of present-day DIC concentrations even at saturating light levels. The P max of the CCM species X. gladiata remained constant across water motion and pH treatment conditions, even when those conditions would theoretically reduce the availability of DIC at the alga's surface. This suggests that X. gladiata does not significantly down-regulate the activity of its CCM when the availability of DIC decreases, instead maintaining a high photosynthetic rate under a wide range of conditions. In contrast, the non-CCM species, Hymenena palmata did show evidence of carbon limitation in terms of its P max . Under present-day pH conditions (pH 8.1), H. palmata (non-CCM) exhibited a lower P max within the low-mixed compared to the high-mixed treatment, a finding consistent with Kübler and Dudgeon 9 who showed with a model that thick boundary layers can limit net photosynthetic rates of non-CCM macroalgae. Interestingly, when the proportion of CO 2 in the total DIC was increased in the low-pH treatment, the P max of H. palmata within the low-mixed treatment matched that of the high-mixed treatment. Unlike X. gladiata, it appears that the photosynthetic rate of this non-CCM red species depends upon the availability of CO 2 , indicating a greater sensitivity to external CO 2 availability than the brown macroalgae with an active CCM.
It was hypothesised that regulation of CO 2 uptake by water motion would be counteracted by a lower pH that increases the concentration of CO 2 . This positive response to increasing CO 2 availability, however, is dependent on the species sensitivity to the external CO 2 concentration. This sensitivity seems to be related to the absence of a CCM and the local environment that the organism inhabits (i.e., habitats with low mixing). The low pH treatment did, however, have an alternative negative effect on the CCM species, X. gladiata, with this species exhibiting a decrease in DIC uptake efficiency in the lower pH treatment compared to the ambient pH, both in www.nature.com/scientificreports/ high-and low-mixed conditions. Johnston and Raven 52 observed a similar finding when Fucus serratus (brown) showed a reduced affinity of CO 2 when grown in high CO 2 conditions 52 . This could be attributed to the CCM or another physiological process of the algae being repressed, however, it is unclear whether this is caused by the low pH or the high CO 2 32,52 . Indeed, a low pH means that there is an increase in the [H + ] within the seawater. As H + plays a vital role in the regulation of cellular homeostasis, an increase in [H + ] could impact metabolic processes and CCM activity of certain species of macroalgae 32 . However, not all CCM macroalgal species show this response and other species with CCMs have shown significantly positive 32 or no photosynthetic or growth response 46 to increasing pCO 2 and lowering pH.
Water motion influences several factors in the coastal environment, impacting not only solute exchange but also the sediment dynamics, the light environment, temperature and grazing 14,25 . Coralline algae were significantly more abundant at the wave-exposed sites compared with the wave-sheltered sites. This is most likely due to more fine sediment and turfing algae being present at the sheltered sites, where they were observed smothering the low-lying encrusting coralline algae. Coralline algae typically possess CCMs and will, to varying degrees, also assimilate CO 2 via diffusion 40 , so are not expected to be limited by DIC. Bergstrom et al. 40 investigated the carbon uptake response of crustose coralline algae to OA and found that only one out of six coralline algae species down-regulated its CCM under elevated pCO 2 , suggesting that most coralline species will not have an obvious positive benefit from increasing CO 2 availability.
Given the contrasting responses to changing CO 2 availability in both the field observations and the laboratory experiments, we provide evidence that macroalgal species have differing levels of sensitivity to the external CO 2 availability. In general, browns and green macroalgae with CCMs are resilient to changing CO 2 availability and maintain an active CCM in a wide range of conditions, while red fleshy algae that show non-CCM behaviour are more sensitive. This sensitivity indicates that within present-day low-mixed environments, large diffusive boundary layers could impede the delivery of CO 2 to the surface of non-CCM algae. We suggest that the distribution of non-CCM species in shallow coastal areas is in part regulated by water motion. However, this probably only has a minor effect on the overall structure of shallow macroalgal assemblages in present day CO 2 conditions in southern New Zealand, as the majority of algae in both the wave-exposed and wave-sheltered sites possess values of δ 13 C indicative of an active CCM 20 . In deeper, low-light field sites that are dominated (> 80%) by non-CCM species (e.g. southeastern Tasmania 34 ), however, water motion is likely to strongly regulate the distribution of macroalgal species. Macroalgal species that currently experience carbon limitation, and which are sensitive to external CO 2 availability for carbon acquisition, could benefit from predicted increases in seawater concentration of CO 2 caused by the increase of atmospheric CO 2 from anthropogenic emissions 53 . Although since the majority of macroalgal species possess a CCM 11,41 and do not appear to be limited under the present-day availability of DIC, a significant change in community structure is unlikely.

Methods
Site information. Surveys of the macroalgal assemblages and tissue collections were conducted within kelp forest habitats between 3 and 5 m below mean low water on the south-eastern coast of the South Island of New Zealand (45°45′27"S, 170°41′37"E). Four wave-sheltered sites within Otago Harbour and three wave-exposed sites along the open coast just north of the harbour were selected (Fig. 5). All sites are linked by the same bulk water, and there are no significant freshwater outflows within the area that might influence one site differently. Macrocystis pyrifera (Ochropyta, Phaeophyceae, Laminariales) was the dominant kelp canopy at all sites. The presence of the bull kelp Durvillaea spp. (Ochrophyta, Phaeophyceae, Fucales) was used as a biological indicator of wave-exposure at the sites, as this species is exclusive to wave-exposed shores 24 .
The rate of dissolution of gypsum blocks was employed to obtain a proxy for the rate of mass transfer of solutes at the study sites 54,55 . Dissolution of gypsum blocks is mediated by the thickness of boundary layers and is driven by advection processes, thereby providing a useful assessment of relative differences in mass transfer rates of solutes between the different sites. The gypsum blocks were made following Gerard and Mann 55 and were deployed by attaching the blocks to an anchored float, where they were positioned 0.5-1 m below the water surface and on the outer edge of the reef at each site. The blocks were deployed on three separate occasions during April 2010 (autumn) for 48-54 h. The hourly dissolution rate was calculated from the before and after dry weight. It was ensured that more than 30% of the original weight of the gypsum blocks remained, so that the rate of dissolution could be considered linear.
Light attenuation, temperature, salinity and nutrient concentrations were measured on three separate days at each of the sites between April and May 2010. To obtain an estimate of the water turbidity at the sites, light attenuation was measured with a LI-COR LI-193 Spherical Quantum Sensor. Single light measurements were taken every meter and the exponential rate of decay of light with depth was calculated as a slope of the linear regression of the natural logarithm on depth 58 . The temperature and light measurements were compared with the extensive light dataset of the area conducted by Pritchard et al. (2013) to ensure the values were representative of the area. Seawater salinity was determined with a Sper Scientific © Salinity Refractometer. Triplicate water samples were collected from the surface at each of the sites in April 2010 and were filtered through Whatman™ GF/C filters. The samples were frozen before the concentrations of nitrate, ammonium and phosphate were analysed on a Lachat QuickChem © 8000 automated ion analyser using standard methods 59 within a month of collection.

Surveys of macroalgal assemblages and δ 13 C analysis.
Benthic surveys were carried out at each of the seven sites during February (summer) 2010 and were completed within one month. 20 A4-sized (210 × 297 mm) underwater photo-quadrats were taken at each of the sites, using a Canon Powershot G10 digital camera in an underwater casing. The photo-quadrats were taken from random positions along a 30 m transect that was haphazardly placed parallel to the shore at a depth of 1-3 m. Percent cover of species and functional groups was Scientific Reports | (2022) 12:21947 | https://doi.org/10.1038/s41598-022-26517-z www.nature.com/scientificreports/ calculated with the random point count method, with 50 random points overlaid on each of the photo-quadrats using Coral Point Count 60 . Macroalgae were grouped into their phylum or class i.e. Chlorophyta (greens), Rhodophyta (reds), or Ochrophyta, Phaeophyceae (browns). The browns were further split into filamentous or leathery, and the reds split into fleshy and corallines, to take into account the different ecological roles these groups play within a community. Corresponding tissue samples were taken from the same location on three replicate individuals of each of the macroalgal species observed at the seven sites were collected at the same time as the community surveys. The thallus samples were cleaned of epibionts, rinsed with filtered seawater, and were identified to species level based on their morphology. The tissue was dried to a constant weight at 80 °C and then ground to a fine powder. δ 13 C (‰) was determined with the ground samples using a CE NA1500 Elemental Analyzer (Carlo-Erba instruments) interfaced to a Europa Scientific 20-20 update continuous flow mass spectrometer (Department of Chemistry, University of Otago, New Zealand). Corrections for drift were made automatically every 5 samples from a standard ethylenediamine tetraacetic (EDTA) with a known isotopic ratio. The δ 13 C values were expressed as ‰ and represent the relative 13 C: 12 C of the sample against a standard ( 13 C: 12 C of Pee Dee Belemnite) multiplied by 1000. The macroalgal species present in each of the assemblages surveyed were grouped into one of three carbon-use strategy groups as used by Hepburn et al. 20 : CCM (δ 13 C ≥ − 29‰), non-CCM (δ 13 C < − 29‰) and coralline algae. There are a number of other potential reasons for low δ 13 C in marine macroalgae, including refixation of respired CO 2 , a 10.72‰ more negative δ 13 C of CO 2 produced from HCO 3 − during intracellular conversion by carbonic anhydrase and exposure to a more negative δ 13 C of external inorganic carbon due to estuarine input 11 . As discussed in Raven et al. 7,11 , these first two cases are unlikely in Rhodophyta macroalgae (the most common non-CCM group in Southern New Zealand), while differences in the δ 13 C of external inorganic carbon between the study sites is improbable due to the sites all being linked by the same bulk water with no estuarine or freshwater outputs in the vicinity.
DIC uptake efficiency and O 2 evolution. In a subsequent laboratory experiment, we examined the effect of water motion and pH on DIC uptake by two species, a CCM and non-CCM, that were dominant in the field surveys. Approximately 1 g thalli samples of Xiphophora gladiata subsp. novae-zelandiae E.L.Rice (brown CCM species) and Hymenena palmata (Harvey) Kylin (red, non-CCM species) were collected from the same location on thalli of different replicate individuals at 1-2 m depth from Butterfly Bay, Karitane, Otago (Fig. 5) between April and May 2011 (Autumn). All thalli samples were taken from the top of the sub-canopy, with only a M. pyrifera canopy above them. Thalli were kept in seawater in a cool box for transport back to the laboratory www.nature.com/scientificreports/ 30 min away. At the laboratory, the thalli were gently cleaned of epiphytes and sediment, and were acclimated overnight in filtered seawater to facilitate wound healing. The thalli were stored for up to 4 days in filtered seawater at a temperature of 12 °C with constant bubbling provided by an air bubbler stone. TrueLumen™ 460 nm Actinic T5HO lamps on a 12:12 light dark cycle provided 100 µmol photons m −2 s −1 of light to the stored samples. The seawater used in the experiments was collected from the Portobello Marine Laboratory (University of Otago) and filtered through a 0.5 μm pore size Filter Pure ® polypropylene spun melt filter and then sterilised with an Aquastep ® 25 W Ultraviolet Sterilizer. A full factorial design was used to measure the effect of two levels of pH (8.1 and 7.6) and two levels of water motion (high-and low-mixed conditions) on the uptake of DIC for photosynthesis in X. gladiata and H. palmata, with 4 replicates. The IPCC's RCP 8.5 emissions scenario 18,19 with a pH of 7.6 by 2100 was used as an extreme future ocean acidification scenario. Experiments were conducted in 150 mL chambers under a photosynthetically saturating irradiance of 300 µmol photons m −2 s −1 at 12 °C. Light was delivered by a Kodak Carousel ® S-AV2000 slide projector. To maintain the treatment pH throughout the experiments and thus a constant ratio of CO 2 :HCO 3 , Tris buffer was added to UV filtered seawater at a concentration of 26 mM. Tris buffer has been reported to affect the photosynthesis of the macroalgae Saccharina latissima (formerly Laminaria saccharina) (Ochrophyta, Phaeophyceae, Laminariales), by interfering with the method of HCO 3 − utilisation this alga possesses 61 . However, the effect of Tris buffer on the photosynthetic ability of X. gladiata and H. palmata was tested by measuring the O 2 evolution of individuals (n = 3) in filtered seawater with and without TRIS, and no detectable effect was found.
DIC was removed from the Tris-buffered seawater by lowering the pH to below 3 with HCl, bubbling with nitrogen gas for 4 h and then raising the pH to either 8.1 or 7.6 with fresh NaOH solution. A stirrer bar was used to modify the level of mixing in the O 2 evolution chambers, which was turned to either a high or low speed to apply a high mixed treatment representing a fast flow environment where boundary layers would be thin or a low-mixed treatment with minimal mixing and where thicker boundary layers would exist. The level of mixing in the chambers was measured by injecting dye into the chambers and measuring the time for the dye to completely disperse into the seawater. The dye was 5 times faster to disperse in the high mixed treatment, indicating a fivefold increase in the mass transfer rate in the high-mixed relative to the low-mixed treatment.
Sixteen individuals of each species were randomly assigned to each of the four treatment combinations, to give four replicates of each treatment. O 2 concentration of the water was continuously recorded with a fibreoptic Ocean Optics ® FOXY-R sensor probe connected to a USB-2000 spectrophotometer. The O 2 concentration of the DIC free Tris-buffered seawater was lowered to 20-40% before each experimental-run by briefly sparging with nitrogen gas (~ 5 min), to reduce any possible effect of photorespiration by the algae. Oxygen evolution measurements were begun once the level of O 2 in the medium was stable (i.e. no photosynthesis was occurring), indicating that individuals had used up any internal inorganic carbon stores or remaining DIC in the seawater medium; this took around 10 min. Every 10 min, 0.05 mL of 1 M NaHCO 3 − was injected into the chamber until the total DIC concentration in the seawater medium was 1.825 mM, thereafter 0.1-0.3 mL additions were made until the total calculated concentration of DIC was 7 mM. This resulted in 12 injections with a total incubation time of 120 min. The photosynthetic rate was measured during each DIC addition, which was used to calculate the DIC-saturated maximum photosynthetic rate (P max ) and DIC uptake efficiency (k 0.5, the concentration of DIC where photosynthesis is half of P max ), which is further detailed in the statistical methods section. The alkalinity and concentrations of CO 2 , HCO 3 − and CO 3 2− of the treatment seawater medium at each addition of DIC were calculated using the known parameters of total DIC, temperature (12 °C), salinity (35 psu) and pH (7.6 or 8.1). The calculations were done using the constants of Mehrbach et al. (1973) refitted by Dickson and Millero (1987), with the R package 'seacarb' 42 and are presented in Fig. 4.
The oxygen evolution measurements were standardised to the algal fresh weight. A linear regression was applied to these standardised oxygen evolution measurements between each addition of NaHCO 3 − to calculate the photosynthetic rate (µmol O 2 g −1 FW s −1 ) of the individuals at each DIC concentration.
Statistical methods. Differences in functional group percent cover and the abundance of each carbonuse strategy, were tested between each site and the two levels of water motion with a nested ANOVA using the R statistical software platform 62 . Each site was nested within either the wave-sheltered or wave-exposed water motion classification factors. Sites were considered replicates of the two levels of water motion, giving 4 wavesheltered replicates and 3 wave-exposed replicates. The variability between the sites was of interest and so a nested ANOVA test was used to provide a test of comparison between the two water motion classifications, but also among the seven sites taking into account the water motion at the sites. One-way ANOVAs were run for the δ 13 C values of any species occurring at both the wave-exposed and wave-sheltered sites, to examine whether they were significantly different between water motion regimes. A p-value equal to or less than 0.05 was considered significant. No transformations of the data to meet the requirements of these parametric tests were required.
A Michaelis-Menten curve was fitted to plots of photosynthetic rate vs DIC concentration. This model is fitted to data by optimising the photosynthetic maximum (P max, the DIC-saturated maximum photosynthetic rate) and DIC uptake efficiency (k 0.5, the concentration of DIC where photosynthesis is half of P max ). A maximum likelihood, non-linear mixed effects modelling approach, using the R statistical software platform 62 and the lme4 package 63 was used. This approach allowed a single model for each species, whilst accounting for variability in P max and K 0.5 , by including incubation as a random factor. Within species contrasts on P max and k 0.5 of X. gladiata and H. palmata were made between the four treatment combinations (low or high mixed treatment with ambient or low pH) using simultaneous t-tests, using the R package multcomp 64